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T  Reduction 

Jiced  systems,  as  defined  here,  are  finite-dimensional  autonomous 
continuous-time  dynamic  systems  governed  by  equations  of  the  form 

(t)  ■  f(x(t));  xQ(to)  *  x0eR"»  t  1  t0»  where  f  :  Rn  Rn  is  piecewise- 
constant  with  discontinuities  only  on  the  surfaces  where  one  or  more  co- 
ordinates  of  Rn  take  integer  values.  A  diced  system  in  R  is  very  easy  to 
illustrate:  the  plane  can  be  divided  into  a  uniform  grid,  and  within  each 
square  a  vector  representing  the  magnitude  and  direction  of  f  is  shown 
(Figure  1). 


Existence  and  uniqueness  of  a  solution  for  any  fixed  initial  state, 
xq,  can  be  studied  using  a  generalization  of  the  method  introduced  by 
Filippov  [1];  trajectories  may  exhibit  sliding  mode  segments  and  higher- 
order  non  differentiable  behavior  as  illustrated  in  Figure  1.  In  order 
to  obtain  existence  of  solutions, 
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multivalued  extensions  of  f  onto  its  discontinuity  surfaces 

are  required.  Every  trajectory  can  be  represented  by  a  sequence  of 

transition  points  and  times,  {x(t^,t.}. 

Definitions  of  various  types  of  stability  and  instability  can  be  con¬ 
structed  from  an  examination  of  the  invariant  limit  sets  [2]  of  the 
trajectories.  For  diced  systems,  the  range  of  asymptotic  behavior  of 

trajectories  starting  from  different  initial  conditions  can  be  exceedingly 

.  2 
rich;  Figure  2  illustrates  some  patterns  in  R  .  The  possibility  of  approx¬ 
imate  global  stability  analysis  using  nondeterministic  automata  is  examined 
and  its  limitations  are  shown. 

In  practice,  diced  systems  might  be  viewed  as  approximations  of  con¬ 
tinuous  or  discontinuous  systems.  In  the  former  case,  for  instance,  we 
might  seek  the  best  piecewise-constant  (finite-element)  approximation  to  a 
continuous  system.  Wang  [3]  has  presented  an  application  of  this  type  for 
solving  partial  differential  equations.  In  the  latter  case,  a  s  ate  space 
diffeomorphism  might  be  used  first  to  transform  the  discontinuities  of  a 
system  to  lie  along  coordinate  axes,  and  then  a  diced  approximation  could 
be  developed  which  would  preserve  the  discontinuous  behavior  of  such  systems. 

The  potential  practical  advantages  of  diced  approximations  lie  in  a  reduction  of 
information  storage  required  to  characterize  a  system  and  the  possibility  of 
assessing  its  approximate  asymptotic  behavior  without  a  detailed  simulation. 

For  example,  at  the  time  of  a  known  failure  of  a  power  system,  it  is 
often  desirable  to  predict  the  long-term  consequences  of  various  control 
strategies  so  that  an  operator  can  decide  among  them.  Yet  the  system  is  too 
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big  to  store  all  possible  consequences  in  advance.  A  practice  which  has 
thus  been  followed  in  some  cases  [4]  is  to  run  a  simulation  "faster  than 
real-time"  for  each  control  strategy.  While  the  issue  of  approximation 
accuracy  is  not  treated  here,  the  results  suggest  that  significant  economy 
of  real-time  computation  might  be  achieved  by  approximating  the  dynamics 
of  a  diced  system.  However,  they  also  suggest  that  the  patterns  of  stability 
and  instability  exhibited  by  such  discontinuous  systems  may  be  highly  com¬ 
plex  and  that  analytical  methods  are  not  likely  to  yield  clear-cut  predictions 
about  global  stability. 
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II.  Preliminaries ,  Notation 

Let  i  *  [ij, . . . ,in] eZn  be  a  multi-index  on  the  n-tuples  of  integers 
(Z).  Let  b  *  [bj, . . . ,bn] eBn  represent  an  n- tuple  of  binary  numbers 
(B  ■  {0,1}).  Let  X^(x):  Rn  <*■  R  be  the  characteristic  function  of  the  open 
set  {x  »  [x^,...,x  ]eRn  |  ifc  <  x^  <  i^+1,  k  *  l,2,...,n}. 

Definition;  A  diced  initial  value  problem  (DIVP)  is  specified  by  a 
system  of  ordinary  differential  equations 

x{t)  *  f(x(t))  ;  x(tQ)  =  xQeRn  ;  t  >  tQ  (2.1) 

where  f  :  Rn  ■*  Rn  has  the  particular  form 

f(x)  *  l  fbi  X^OO  ;  b  »  0£Bn  (2.2) 

ieZn 

and  e  Rn  for  each  multi- index  i.  * 

The  surfaces  of  discontinuity  of  f  may  be  classified  by  their  dimension. 
Let  1(b)  :  Bn  (l,...,n)  be  a  function  denoting  the  number  of  "l'"s  in  the 
binary  n-tuple  b.  For  fixed  i  s  Zn,  consider  the  sets 

hi  a  {xeR°  I  h  <  *k  <  h*i  if  bi  *  0 

ik  ■  x  if  bk  »  1,  k  »  1,2,..., n}  (2.3) 

These  may  be  viewed  as  the  set  of  submanifolds  "attached  to"  the  point  x*i . * 
For  example  is  the  interior  of  the  n-dimensional  cube  indexed  by  its 
vertex  at  x«i;  S.^  (the  shorthand  1  denoting  b  *  [1,1,...,  1])  is  the  single 
point  x»i.  The  submanifolds  of  dimension  p  associated  with  x*i  are 


•The  obvious  injection  of  the  integers  into  the  reals  is  implied. 


e. 
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S?  *  {S^  |  1(b)  »  n-p}  p  ■  0,1,..., n.  (2.4) 

This  notation  provides  a  compact  classification  of  all  of  the  subsets  of 
Rn  which  are  of  interest. 

In  Section  III,  conditions  for  well-posedness  of  a  OIVP  are  examined. 

This  is  done  by  extending  f  to  its  discontinuity  surfaces 

(from  {foi>,  we  generate  {fbi>.  b  i  0eRn).  Then  a  constructive  procedure 

can  be  used  to  generate  solutions  x(t)  *  $(t,to,xQ)  ^or  each  x0e  R”>  t  e  R 

and  hence  to  define  the  transition  map  4>:  R  x  R  x  Rn  Rn.  Let  X  denote 

the  function  space  in  which  trajectories  are  defined.  This  leads  to  the  following 

Definition:  A  diced  system  is  an  autonomous  dynamical  system  (X,Rn,$) 

(See. [S]). 

Stability  has  been  viewed  as  a  qualitative  property  of  a  dynamical 
system,  and  concerns  the  asymptotic  behaviors  of  trajectories 
x(*)  »  '♦'(•, t  ,XQ)  as  xQeX  is  varied.  Stability  of  diced  systems  is  dis¬ 
cussed  in  Section  IV.  Two  useful  notions  will  be  those  of  the  positive 
limit  set  and  the  invariant  set  f2]> 

Definition:  The  set  (1  C  Rn  is  invariant  with  respect  to  the  system 
x(t)  *  f(x(t),t)  if  for  any  xQ  e  £2  there  is  a  tQ  such  that  the  motion 
5(t.t  ,x  )  belongs  to  8  for  all  t  >  tQ. 

Definition:  The  set  H  C  Rn  is  called  the  positive  limit  set  of  a 

bounded  motion  <Kt;to,xQ)  if,  for  any  point  pen,  there  exists  a  sequence 

of  times  {t  }  tending  to  infinity  as  so  that 

n 

am  I  "  pH  *  0  (2,5) 

In  applying  these  definitions  it  will  be  useful  to  recall  that  a 
function  x(t)  is  periodic  T  >  0  if  x(t)  ■  x(t*T)  for  all  t;  "the"  period 
of  a  periodic  function  is  defined  as  the  least  T  for  which  this  equality  holds. 


\ 
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III.  Existence  and  Uniqueness 

Consider  the  DIVP  (2.1), (2.2).  Defining  solutions  within  the  cubes 
by  integration  is  entirely  straightforward;  all  difficulties  arise  in 
attempting  to  extend  solutions  across  the  discontinuity  surfaces  of  f;  in 
general,  there  is  no  unique  continuation.  Various  possibilities  are 

(a)  To  restrict  the  class  of  f  so  that  continuations  are  always  unique 
(this  is  very  restrictive  indeed,  and  essentially  eliminates  many 
interesting  phenomena  from  consideration) . 

(b)  To  eliminate  the  non-continuable  surfaces  from  the  domain  of  f;  how¬ 
ever,  then  all  points  on  all  trajectories  leading  to  such  surfaces 
must  also  be  eliminated,  and  a  large  part  of  the  original  domain  of 
definition  may  ultimately  be  excluded. 

(c)  To  choose  an  ad  hoc  rule  for  continuation  of  solutions;  however,  it 
proves  difficult  or  impossible  to  do  this  in  a  self-consistent  and 
unbiased  manner. 

A  fourth  alternative  has  been  selected  here: 

(d)  To  sacrifice  uniqueness  and  continue  all  solutions  through  a  dis¬ 
continuity. 

In  this  way  a  viable  deterministic  existence  theory  can  be  developed,  at 
the  cost  of  considering  a  countable  number  of  alternative  solutions.  A 
"physical"  justification  for  adopting  this  approach  is  that  in  the  presence 
of  small  perturbations  of  the  initial  conditions,  a  solution  near  to  at 
least  one  alternative  solution  will  occur. 

A  constructive  procedure  is  given  for  defining  solutions.  To  simplify 


its  presentation,  a  multivalued  continuation  of  f  to  the  surfaces  S^,  b  /  0 
is  first  defined.  Initially,  f  is  specified  on  the  submanifolds  S?  *  {SQ^} 
of  dimension  n.  The  continuation  proceeds  recursively  to  submanifolds, 

S?  of  dimension  n-1,  n-2,...,0.  Recall  that  S?  is  the  point  set 
(x  £  Rn  I  x  *  i^,  i^  an  integer).  Notationally,  a  single  valued  f^  will 
not  be  distinguished  from  a  multivalued  f^,  the  implication  being  that 
the  prescribed  rule  is  applied  to  each  possible«value  of  f^  in  turn,  and 
the  set  of  all  results  is  retained.  Let  p  *  n.  Suppose  f^  are  known  on 
S?,  p  £  q  £  n.  Then  f^  can  be  extended  to  S?"1  as  follows,  for  each  i  z  Zn 

Suppose  £  S?”1 .  Let  indices  j  j. . .  (p-1)  ^enote  t^ie  ordered 

nonzero  positions  of  b,  i.e.,  b.  *  1,  k  3  1, . . . ,n- (p-1)  and  b.  *  0  other- 

Jk  J 

wise.  The  neighborhoods  of  of  dimension  q,  p  £  q  £  n,  can  be  defined 
as  follows.  For  q  »  n,  consider  all  indices  I  formed  by  decrementing 
ij^  by  one  for  any  subset  of  the  subindices  k  *  1,.. .n-(p-l),  including 
the  null-set;  then  Sgj  £  S?  is  a  neighborhood  of  where  S  *  0.  For 

q  ■  n-1,  consider  all  values  b  having  a  single  "one"  in  one  of  the  positions 
j 1 . . . ^n  (p  i)  and  for  each  B,  form  1  from  the  remaining  n-(p-l)-l  indices 
as  above;  then  Sgr  £  is  a  neighborhood  of  S^.  For  q  *  n-2,  con¬ 

sider  all  values  B  having  "ones"  in  any  two  of  the  positions  j^, • • • »^n-(p-l) 
and  from  each  6  form  1  from  the  remaining  n- (p-1) -2  indices  as  above;  then 
Sg>  £  is  a  neighborhood  of  .  This  procedure  is  continued  until 
q  »  p. 

The  values  of  f.  .  on  S.  .  £  S?”*  are  determined  from  the  values  of 
01  01  1 

fgj  on  each  of  its  neighborhoods  Sg*  £  s3,  p  £  q  £  n.  It  is  thus 
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sufficient  to  give  the  procedure  for  determining  f^,  assuming  that  these 

values  on  higher-dimensional  submanifolds  are  known  (i.e.,  the  values  can 

be  determined  recursively) .  Define  Sgr  to  be  an  input  submanifold  to 

if  (fr-r)  *  0  for  all  Z  such  that  5»  »  1,  and  for  all  remaining  Z  in  the 
bl  l  1 

set  (p-i)*  <  0  for  those  *  such  that  H  a  H»  while 

z 

>  0  ^or  ^ose  *  suc^  tbat  H  *  ig-1*  Define  Sgr  to  be  an  output  sub¬ 
manifold  if  *  0  for  all  Z  such  that  *  1,  and  for  all  remaining 

l  in  the  set  jj,.  •  *jn-(p-l) '  ^fbi\  —  0  for  those  *  such  that  H  *  while 
tfbi\  -  0  for  those  A  such  that  H  *  ig-1.  Note  that  those  sets  for 
which  (fgj)  J*  0  when  6^=1  need  not  be  considered.  So  long  as  the  set  of 
output  submanifolds  of  is  non-empty,  f^  is  assigned  the  set  of  all 
values  fgj  on  the  output  submanifolds.  If  the  set  of  output  submanifolds 
is  empty,  is  a  generalized  sliding  surface.  Consider  fg.  on  Sg-  t  s5  in  the 
input  set.  If  this  set  is  empty,  set  *fgr  *  0.  Recall  that  Sgr  e  S?  is 
formed  by  keeping  i  unchanged  in  all  but  one  position,  say  jk,  of  b,  so 

b  *  [b,,...b.  .,  0,  b.  .  ,...b  ]  and  either  i  =  i  or  i  *  [i. , . . . ,i .  -1 , 

Jk  Jk  n  1  Jk 

i.  .,..., i  ].  Thus  there  are  a  maximum  of  2(n-(p-l))  surfaces  in  this 

h  1  n 

subset  of  the  input  set.  These  surfaces  are  considered  in  pairs  to 

determine  the  admissable  values  of  fg^;  using  the  example  above,  if  is 

in  the  input  set  then  (fg^)j  <  0  and  if  Sgr  is  in  the  input  set  (fg?) ^  >  0. 

k  k 

If  both  elements  are  members  of  the  input  set  then 

fM  *  [(fSlVSi  -  -  CfjiJjJ  (3.1, 

while  if  only  one  is  in  the  input  set,  let 
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The  set  of  possible  values  of  f^  on  a  generalized  sliding  mode  is  coupleted 
by  considering  each  Sg»  £  sE  in  this  manner.  In  all  such  cases,  ^i,; 
k  »  1, . . . ,n-(p-l)  are  zero,  so  that  further  motion  occurs  on  itself. 

Thus,  the  procedure  for  extending  the  function  f  to  all  of  Rn  is  com¬ 
pleted.  The  complexity  of  the  procedure  arises  from  the  large  number  of 
possibilities  which  can  arise.  A  number  of  such  special  cases  are  illus¬ 
trated  on  Figure  3.  Evidently,  the  procedure  for  extending  f  is  not  the  only 
one  which  could  be  devised.  In  the  next  step,  construction  of  solutions, 
however,  it  will  become  apparent  that  the  underlying  principle  has  been 
to  define  f  in  a  manner  which  preserves  all  trajectories  that  might  arise 
from  each  initial  condition. 

Let  xq  £  Rn  be  given  as  the  initial  condition  of  (2.1)  at  t  =  tQ;  let 
£  S?  be  the  smallest  submanifold  containing  xq.  Let  f^  denote  one  of 
the  extended  values  of  f  on  S^.  Define 


4>(t,t  ,x  )  *  x  +  f.  .  (T-t  )  ;  t  <  t  <  t 

oo  obio  o  —  l 


(3.2) 


The  time  t^  is  defined  as  follows:  for  each  i  such  that  (f^j^  *s  nonzer°. 

let  (tj)^  denote  the  first  T  >  t  such  that  ^  is  an  integer; 

then  tj  *  min[(t1)^]  and  Xj  *  <Kti*t0»x0i*  If  then  *  «  and 

x  *  x  ,  and  this  soLution  terminates.  Otherwise,  x.  defines  new  values 
lo  l 

of  b,i,  and  p,  and  the  solution  process  continues: 

♦<T- W  *  *k  *  fbi(T-V  ;  *1 « T  i  Vi  (3-3) 


I 
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On  those  surfaces  where  f^  is  multivalued,  each  possibility  must  be 
examined  in  turn;  in  this  sense,  is  also  multivalued.  Each  trajectory 
pieced- together  in  this  fashion  can  be  summarized  by  a  sequence 
{.xk,tk},  k  =  0,1,...  in  some  cases,  these  sequences  are  finite  and  in 
other  cases  infinite.  By  inspection  of  {x^}  alone,  a  corresponding 
sequence  of  regions  where  is  the  minimal  submanifold  con¬ 


taining  x^,  can  be  constructed. 

A  solution  of  (2.1),  (2.2)  is  then  defined  in  the  obvious  manner,  as 
any  $(t,to,xo)  constructed  by  the  continuation  procedure  (3.3).  It  has 
the  property  that  for  any  finite  admissablek,  $(t,to,xo)  is  piecewise 
continuous  on  [t  ,t,].  This  solution  by  continuation  is  said  to  be 

0  K 

asymptotic  if  lim  tv  *  °°.  An  asymptotic  solution  is  piecewise  continuous. 

- -  k  K 

For  purposes  of  the  present  work,  a  solution  will  be  said  to  exist  if  the 
state-space  continuation  is  asymptotic.*  Asymptotic  solutions  need  not 
be  unique,  but  the  rate  of  growth  in  the  number  of  solutions  can  be  bounded 
as  a  function  of  k,  since  the  maximum  number  of  output  submanifolds  can  be 
bounded  above  for  any  S^.  If  there  is  only  one  asymptotic  solution  through 
(xo,tQ) ,  it  is  said  to  be  unique.  Continuous  dependence  of  4>(t,tQ,xo)  with 
respect  to  xq,  of  course,  is  not  to  be  expected  for  t  >  t^. 


*More  over  if  lim  t.  t* 08 ,  solutions  by  time-continuation  could  be 

defined;  however,  their  properties  will  not  be  explored  here. 


.1  i 


vr 
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IV.  Stability 

The  usual  definitions  of  stability  presuppose  a  solution  which  is  well- 
posed  in  the  sense  of  existence,  uniqueness,  and  continuous  dependence  on 
the  initial  data.  Diced  systems,  in  general,  do  not  posses  the  last  two 
properties.  One  alternative  is  to  nevertheless  use  the  standard  notions 
of  stability,  restricting  their  domain  of  application  to  those  initial 
states  for  which  the  usual  notions  of  well-posedness  are  (locally)  satisfied. 
Unfortunately,  the  set  of  such  initial  states  appears  quite  difficult  to 
characterize  and  thus  imposes  an  awkward  restriction  on  the  applicability 
of  this  alternative. 

Another  alternative,  introduced  here,  does  not  impose  such  restrictions, 

but  weakens  the  notion  of  stability  that  is  employed.  Stability  is  viewed 

as  a  qualitative  property  of  a  trajectory,  and  a  system  is  then  said  to  be 

stable  when  all  of  its  trajectors  share  this  property. 

Definition:  The  motion  of  a  diced  system  (2.1),  (2.2)  initiated  at 

(t  ,x  )  is 
o’  o 

M(tQ,xo)  *  {$(t,tQ,xo),  t  £  t  |  p  is  a  transition  function 

initiated  at  (t  ,x  )} 
o  o 

which  is  the  set  of  all  trajectories  originating  at  (tQ,xo). 

Definition:  The  motion  M(tQ,xo)  of  a  diced  system  is  said  to  be 
(a)  Bounded  in  magnitude  if  there  is  a  constant  <5  >  0  such  that 

max  {sup  I  |<|>(t,t  ,x  )|  |}  <  P 

^M(to.xo)  t>tQ 
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(b)  Bounded  in  cardinality  if  there  exists  a  constant  N  such  that 

sup  {cardinality  of  $(t,t  ,x  ) }  <  N 

t>t  00 

—  o 

The  concepts  of  boundedness  in  magnitude  and  cardinality  are 
independent.  In  both  cases,  the  only  difficulties  occur  at  t  ■+■•,  since 
(a)  any  $(t,to,xQ)  is  by  construction  bounded  for  all  finite  t,  and  (b) 
the  cardinality  of  <fr(t,to,xo)  is  finite,  by  construction,  for  all  finite 
t.  The  following  propositions  are  almost  immediate. 

Proposition  1:  In  (2.1),  suppose  ||fQ^||  <  F  for  all  i,  then 
|  |<fr(t,t0,x0)  -  xQ  |  I  <  F(t-tQ)  for  all  <j>  e  M(tQ,xo) . 

Proof :  The  extension  of  f _ .  to  f,  .  always  guaranteed  that 

-  oi  bi 

I  ifbiM  <  F«  and  the  construction  procedure  (3.3)  guaranteed  that  the 
estimate  of  the  proposition  held  for  each  t.  q,e.d. 

Proposition  2:  Let  |i|  *  |i^|+. . .+|in| .  Suppose  for  system  (2.1) 

there  exists  B  >  0  such  that  for  all  |i|  >  B,  and  k  »  l,...,n,  (foi)k  <  0. 

Then  M(t  ,x  )  is  bounded  in  magnitude, 
o  o 

Proof :  For  any  i  such  that  |i|  >  B,  every  set  &.  contains  output 
submanifolds  with  the  same  |i|  or  smaller  |i|,  and  input  submanifolds  with 
the  same  |i|  or  larger  |i|;  furthermore,  always  outputs  to  Sq_  with 
| l |  <  | i | .  Thus  the  construction  process  cannot  terminate  for  | i |  >  B ,  and 
for  such  i,  |i|  is  reduced  at  least  once  every  n  intervals;  hence  every 
solution  satisfies  |$(t,tQ,xo)|  <  B  for  t  sufficiently  large.  Thus 
M(tQ,xQ)  is  magnitude-bounded. 


i 


-15- 


Proposition  3:  Suppose  that  for  every  i  e  Zn,  b  e  Bn,  has  at  most 
one  output  submanifold.  Then  the  motion  M(tQ,xo)  of  (2.1),  (2. 2)  is  bounded 
in  cardinality. 

Proof:  The  extension  procedure  of  Section  III  shows  that  in  this  case 
f^  takes  the  value  on  its  output  submanifold  or  the  value  zero.  If  a 
trajectory  enters  S^,  it  either  continues  uniquely  to  the  output  sub¬ 
manifold,  or  terminates  at  S^.  In  either  case,  the  cardinality  of  the 
solution  cannot  increase  during  its  construction. 

Thus  there  are  two  notions  of  instability  for  diced  systems:  solutions 
may  become  unbounded  in  magnitude,  and/or  they  may  become  unbounded  in 
cardinality.  This  second  form  of  instability  is  new:  a  trajectory  can 
fracture  and  a  chain  reaction  of  subsequent  fractures  may  ensure — the 
complexity  of  the  process  grows  without  bound. 

Next,  a  notion  of  stability  is  put  forth.  Suppose  that  the  motion 
M(tQ,xo)  of  a  diced  system  is  bounded  in  magnitude  and  cardinality  (or 
simply  ,rbounded") .  Then  a  set  SCR11  consisting  of  a  finite  union  of  the 
submanifolds  is  termed  a  positive  limit  set  of  a  (bounded)  trajectory 
<Kt,t  ,x  )  if  for  any  point  x  t  S,  there  exists  a  sequence  of  times  {t.  }, 
tending  to  infinity  as  k -*■<*>,  so  that 


lim)<KT  t  x  )  -  x(  *  0 

k-w  k  o  o 


(4.1) 


where  )  (  denotes  the  set -membership  metric,  i.e.,  if  x  e  S.^, 


)y-x( 


y  e  S, 

y  i  s, 


bi 


bi 
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In  applying  this  definition,  it  is  important  to  recall  the  standing 
assumption  from  Section  III,  that  all  trajectories  are  asymptotic,  so 
that  such  sequences  {t^}  exist. 

Definition:  A  bounded  motion  M (t  ,x  )  of  a  diced  system  is  termed 

o  o 

pointwise  stable  if  all  trajectories  <{»(t,to,xo)  e  M(tQ,xo)  have  the  same 
positive  limit  set.  The  motion  is  locally  stable  for  xq  e  if  all 
trajectories  $(t,to,x)e  M(tQ,x),  x  e  S^,  have  the  same  positive  limit  set. 
The  motion  is  globally  stable  if  all  trajectories  4>(t,to,x)  have  the  same 
positive  limit  set. 

Concepts  of  uniform  stability  will  not  be  discussed  since  only  time- 
invariant  diced  systems  are  considered  in  the  present  account.*  In  fact, 
the  evaluation  of  stability,  according  to  the  definitions  given,  can  be 
based  merely  on  knowledge  of  the  sequence  {q^}  of  submanifolds  containing 
{x^},  since  it  is  known  from  the  construction  procedure  that  t^+^  >  t^  and 
from  the  asymptotic  assumption  that  lim  t.  »  » .  This  suggests  that  a  wav  to 

k-xa  K 

generate  the  sequence  {o^}  autonomously,  without  explicit  integration  and 
generation  of  {x^,t^}  would  be  particularly  valuable  in  the  assessment  of 
stability.  This  has  not  been  achieved  yet. 

Knowledge  of  the  time- structure  {t. }  of  individual  solutions  can  be 
of  further  value  in  refining  stability  notions.  To  simplify  the  remaining 
concepts  it  is  now  assumed  that  the  trajectories  are  uniquely-defined 


"The  results  could  be  extended  in  this  direction  for  systems  with  continuous 
time-variation;  however  discontinuous ly  time-varying  systems  may  not  be  con- 
timuable,  as  Filippov  pointed  out. 
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(e.g.,  as  occurs  in  Proposition  3)  and  bounded.  Suppose  H  is  a  positive 
limit  set  of  such  a  solution  in  the  conventional  sense  of  Section  II 
(eq.  (2.5)).  Then  in  the  usual  manner  it  can  be  shown  that  H  is  bounded, 
closed,  non-empty  and  invariant,  the  last  property  being  a  consequence  of 
time-invariance.  In  fact,  as  a  consequence  of  finite-dimensionality  of 
Rn,  all  such  solutions  are  asymptotically  almost-periodic  [g] .  Two  cases 
of  special  interest  are  the  asymptotically  constant  (equilibrium)  solution 
and  the  asymptotically  periodic  solution.  These  can  be  identified 
directly  from  the  sequence  {x^,t^}  characterizing  (f>  (t ,  t^ ,  xq)  . 

,t^}  is  finite  of  length  N,  the 
positive- invariant  limit  set  consists  of  one  point,  the  last  value  x^ 

.t^}  is  jointly  periodic  of  period 
m  for  k  >  N,  then  the  positive- inariant  limit  set  is  a  cycle  (closed  curve) 
in  Rn. 

Proof:  For  the  first  case,  note  that  the  construction  procedure 
automatically  defines  t^  *  when  the  sequence  is  finite,  and  this  implies 
a  constant  solution  for  t  t^.  In  the  second  case,  note  that  since 
{x^.t^}  completely  specify  <Kt,tQ,xo),  $  must  be  periodic  of  period 
tk+o’tk*  k  >  whenever 
linear  interpolation  between  these  points) . 

It  is  interesting  to  note  that  for  diced  systems,  the  establishment 
of  an  equilibrium  or  periodic  solution  after  a  finite  time  (t^)  is  often  to 
be  expected  (whereas  this  would  be  considered  exceptional  in  the  case  of 
continuous  differential  equations) ;  however,  in  some  cases  almost  periodic 
solutions  may  also  exist. 


tx^.t^}  is  periodic  (in  fact,  the  solution  is  a 


(for  which  t^  *  °°) .  If  the  sequence  (x^ 


Proposition  4:  If  the  sequence  (x^ 
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V.  Discussion  and  Conclusions 

Tho  present  account  of  the  stability  of  diced  systems  leaves  a  number 
of  original  question  unanswered  and  raises  some  new  ones.  A  study  of  methods 
for  temporal  continuation  of  non- asymptotic  solutions  is  needed;  such 
solutions  may  represent  a  new  sort  of  sliding  mode  which  can  arise  in 
higher  dimensional  spaces,  as  suggested  by  an  example  of  Utkin  [7].  The 
possibility  of  extending  the  techniques  developed  here  to  time-varying 
systems  has  been  mentioned;  Filippov's  general  existence  results  apply  to 
this  problem.  A  study  of  the  partitioning  of  initial  states  which  is  implied 
by  the  proposed  stability  definition  would  also  be  fruitful;  what  properties 
are  shared  by  initial  state  sets  giv  ing  rise  to  the  same  asymptotic  solution? 
In  general,  it  would  appear  that  the  initial  states  within  a  given  region 
can  ultimately  end  up  widely  dispersed.  The  possibility  of  using  an 
automaton  to  simplify  the  propogation  of  solutions  has  also  been  raised. 

The  approximation  of  continuous  systems  by  diced  systems  has  not  been  ex¬ 
plored,  but  under  appropriate  conditions,  a  bound  on  the  approximation  error 
should  be  achievable. 

In  spite  of  the  questions  that  are  unanswered,  some  modest  progress 
has  been  made  toward  defining  the  stability  properties  of  diced  systems. 

First,  a  constructive  continuation  procedure  for  higher  dimensions  has  been 
found;  the  problem  readily  evades  one's  intuition  about  n  ■  1,2  and  even  3 
as  endless  combinations  of  difficult  situations  may  occur.  Second,  a 
compromise  on  the  issue  of  uniqueness  has  been  put  forth:  the  number  of 
admissable  solutions  at  any  finite  time  is  bounded.  Third,  the  concepts  of 
stability  have  been  generalized  to  provide  meaningful  criteria  for  discon- 
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tinuous  systems  of  diced  type. 

Returning  to  the  electric  power  system  example  cited  in  the  opening 
section,  it  would  appear  that  the  implications  of  the  research  might  be 
very  disturbing,  for  two  primary  reasons.  First,  a  new  type  of  instability-- 
an  unbounded  growth  in  the  number  of  possible  solutions  with  time--has  been 
identified.  Second,  and  independently,  the  partitioning  of  the  initial 
state — at  least  in  worst-case  situations- -based  on  asymptotic  properties, 
appears  to  be  very  fine  and  irregular;  thus  a  small  perturbation  in  the 
initial  state  may  give  rise  to  completely  different  asymptotic  behavior 
than  is  found  for  the  unperturbed  initial  state.  Both  of  these  phenomena 
imply  that  the  future  behavior  of  a  diced  system  with  a  (approximately) 
specified  initial  state  may  be  fundamentally  unpredictable;  if  the  long¬ 
term  future  consequences  of  a  present  control  policy  are  unpredictable, 
the  problem  of  choosing  the  best  policy  becomes  more  difficult  and  planning 
must  be  done  with  a  shorter  horizon. 


-20- 


References 


[1]  Filippov,  A.F.,  [1964],  "Differential  Equations  with  Discontinuous 
Right-Hand  Sides,"  Amer.  Math.  Soc.  Trans..  Ser.  2,  Vol.  42,  pp. 
199-231. 

[2]  Willems,  J.L.  [1970],  Stability  Theory  of  Dynamical  Systems,  Nelson 
and  Sons,  Ltd.,  London. 

[3]  Wang,  P.K.C.,  [1968],  "A  Method  for  Approximating  Dynamical  Pro¬ 
cess  by  Finite-State  Systems,"  Int.  J.  Control,  Vol.  8,  No.  3, 
pp.  285-296. 

[4]  Ewart,  D.N.,  "Whys  and  Wherefores  of  Power  System  Blackouts," 

IEEE  Spectrum,  April  1968,  pp.  36-41. 

[5]  Willems,  J.C.,  and  S.K.  Mitter  [1971],  "Controllability,  Observability, 
Pole  Allocation  and  State  Reconstruction,"  IEEE  Trans.  Auto.  Control, 
Vol.  AC-16,  No.  6,  pp.  582-S9S. 

[6]  Dafermos,  C.M.  [1974],  "Semiflows  Associated  with  Compact  and  Uniform 
Processes,  Math  Sys.  Theory,  Vol.  8,  No.  2. 

[7]  Utkin,  V.I.  [1978],  Sliding  Modes  and  Their  Application  in  Variable 
Structure  Systems,  Mir  Publishers,  Moscow.  “ 


INC LASS IF IKD 


SECURITY  CLASSIFICATION  O*  This  PAGE  'When  r>mtm  hnf etrJl 


'  REPORT  DOCUMENTATION  PAGE 


RF  AD  INSTRUCTIONS 
HKl'OKK  COMPLFTIM  •  K)RM 


i  A  ‘  •  mu  \  \2  GOVT  accession  noIT  recipient- s  c  at  a"^ 

AfefegTiu  8^-0547,!  i to-AcUifk  a  .  \  -  • 


C  < and  Submf*} 


STABILITY  OF  DICLD  SYSTKMS  0 


L . 


I 


Ag  THQWf  *i 

T.  Ly/johnson 

IV 


ERFORMING  ORGANIZATION  name  and  address 

MIT  f, 

Laboratory  for  Information  £  Decision  Systems  4 
Cambridge , MA  02139 


S  TYPE  of  REPORT  &  PERIOD  COVEREO 

/ 


Interim, 


Dumber 


•  COnTRACi’^  ORRlfT'  NUMBER;  sj 

/s.  '/tn^n- .  -  • 

F4 9620-80-0-0002  Vt 

To  plo«wRMTvfuEjV*M<!)jggf^WWr 


area  a  work  unit  numbers 


61102F  2304 


U  REPORT  D A*nT* 

20 


y->i  / 


11  CONTROLLING  OFFICE  NAME  AND  AOORESS 

Air  Force  Office  of  Scientific  Research/ NM 
Bolling  AFB  Washington ,  DC  20332 


T3  MONITORING  AGENCY  NAME  a  a  DOR  ESSO  t  dttierent  from  Controlling  OUtc0) 


IS.  SECURITY  CLASS  (ot  thi»  report) 


UNCLASSlif  1LD 


iSa.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


1ft  DISTRIBUTION  STATEMENT  <ol  thi  •  Report) 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  ( ot  the  ebmtrmct  entered  in  Block  20,  It  different  from  Report) 


«•  SUPPLEMENTARY  notes 


IB.  KEY  WOROS  (Continue  on  reveree  etde  If  neceeeery  end  identity  by  block  number) 


10  ABSTRACT  (Continue  on  reeeree  etde  If  noeoooory  end  Identify  by  block  number) 

Diced  systems  are  defined  as  autonomous  systems  governed  by  ordinary  differen¬ 
tial  equations  having  discontinuities  (in  Rn)  on  submanifolds  where  one  or  more 
of  the  state  variables  takes  an  integer  value.  Such  systems  may  be  regarded  as 
approximations  of  continuous  systems  or  at  representative  models  of  a  class  of 
discontinuous  systems.  Trajectories  of  such  systems  (for  a  given  initial  state)! 
are  readily  calculated  and  m*y  exhibit  complex  sliding-mode  segments.  Asymptotip 
properties  of  such  trajectories  are  discussed  and  classified.  Motivation  is 
^{iven^iii^term^^£j>b*ervedjpropert^ 


VT5  °iS0  fri 


00  ' 


f  IAN  T» 


un  seine*  tf  »  *ev  u  it  bmoliti 


MNTT  A^TFTFn 


END 


DATE 

FILMED 


DTIC 


